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Site and bond percolation thresholds are calculated for the face centered cubic, body 
centered cubic and diamond lattices in 4, 5, and 6 dimensions. The results are used to 
study the behaviour of percolation thresholds as a functions of dimension. It is shown 
that the predictions from a recently proposed invariant for percolation thresholds are 
not satisfactory for these lattices. 
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1. Introduction 

Percolation problems have a wide range of applicability, and have therefore attracted 
a fair bit of attention over many yearsEl. Nevertheless the percolation thresholds, 
which are among the basic quantities for percolation on lattices, have been calcu- 
lated exactly for only a few two-dimensional lattices. For many other lattices these 
thresholds have been calculated numerically. These numerical values can then be 
analysed, to determine regular behaviour, or trends, as a function of the lattice 
coordination number, dimensionality, etc. A limited review of such efforts has been 
given elsewhereB. 

Recently, Galam and Maugeru proposed an invariant for site and bond percola- 
tion thresholds, p cs and p c b respectively. The proposed invariant reads 

(p cs ) l/a ^ c6 r l/a& = (i) 

where d is the dimension of the lattice, and a s , and 8 are positive constants. 
The lattices studied by Galam and Mauger were divided into two classes, and for 
each of these classes the values of the parameters were fitted. The values for the 
first class were {a s = 0.3670; a& = 0.6897; 8 = 1.3638}, while for the second class 
{a s = 0.6068; a b = 0.9346; 8 = 1.9340}. For the lattices used by Galam and Mauger, 
the numerical results for this invariant are indeed constant within 5%. Although 
the deviations up to 5% cannot be explained by the inaccuracy in the values for the 
percolation thresholds, it is an interesting observation that the above combination 
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Table 1. The site percolation thresholds of ci-dimensional fee lattices as a function of linear lattice 
size L. The values in the last row are results of a fit to the scaling relation \p c (L) — p c (oo)| oc L~ V". 
Error estimates concerning the last digit are indicated between brackets. 

d=3 d=A d=5 d=6 



L 


Pes 


L 


Pes 


L 


Pes 


L 


Pes 


16 


0.2097(2) 


16 


0.0856(2) 


8 


0.0457(2) 


8 


0.0257(2) 


32 


0.2035(2) 


21 


0.0847(2) 


12 


0.0435(2) 


10 


0.0258(2) 


64 


0.2016(2) 


32 


0.0845(2) 


16 


0.0432(2) 


12 


0.0252(2) 


128 


0.2001(2) 


18 


0.0844(2) 


21 


0.0431(2) 


14 


0.0252(2) 


250 


0.1998(2) 


64 


0.0843(2) 


32 


0.0432(2) 


16 


0.0252(2) 


oo 


0.1994(2) 


oo 


0.0842(3) 


CO 


0.0431(3) 


oo 


0.0252(5) 



of percolation thresholds yields almost constant values. Especially the absence of 
coordination number q in the invariant makes it valuable. 

However, the invariant can be tested on more lattices for which the percolation 
thresholds are known. We list data for many lattices in this article, obtained from 
various sources in literature. However, for some of the higher dimensional lattices, 
only site- or only bond percolation thresholds were known. Therefore the 'missing' 
values are computed here, so that the invariant can be calculated for these lattices 
as well. Moreover, it is interesting to study the scaling of percolation thresholds as 
a function of dimension for some of the important lattices, like the face centered 
cubic lattice (fee), the body centered cubic lattice (bec) and the diamond lattice. 

This paper is organised as follows. In section 2, definitions are given for the 
fee, bec, and diamond lattices in higher dimensions. The cluster algorithm and 
some implementation issues are discussed in section 3, after which the results are 
presented in section 4. Finally, section 5 contains a discussion of the results. 

2. Lattice definitions in higher dimensions 

For the definitions of lattices in d dimensions, we refer to Conway and SloancB fox 
the face centered cubic and body centered cubic lattice, and to Van der MarckB 
for the diamond and Kagome lattice. For completeness, we briefly describe here 
the generalisation of the fee, bec and diamond lattices to higher dimensions. First 
we need to have, for each lattice, a set of d independent lattice vectors. A point 
is called a lattice site if and only if it is an integer combination of these vectors. 
Secondly, we want to know the neighbours of a site, described in terms of the lattice 
vectors. We will denote an orthonormal basis of d-dimensional space by Xj. 

The d-dimensional fee lattice is the set of points in Z5 n , for which the sum 
the coordinates is even. Conway and Sloane use the notation 'D^ for this latticel 
(Note that J} d is the 'closest packed' lattice for d — 3, 4 and 5, but not for higher 
dimensions!! ) Each site has 2d(d — 1) neighbours at a relative location ±Xj ± Xj, 
for i,j = 1, . . . ,d (i ^ j). As a set of lattice vectors one can choose fj = Xi + x«. 
The neighbours of a site are given in terms of these lattice vectors as 



±(f*-fj-) 
±(f 4 + f, - fi) 



foralH,j = l,...,d (i^j). (2) 
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Table 2. The site percolation thresholds of d-dimensional bcc lattices as a function of linear lattice 
size L. The values in the last row are results of a fit to the scaling relation \p c (L) — p c (co)| oc L~ V". 
Error estimates concerning the last digit are indicated between brackets. 





d = 3 




d = 4 




d=5 




d = 6 


L 


Pes 


L 


Pes 


L 


Pes 


L 


Pes 


16 


0.2593(2) 


16 


0.1078(2) 


8 


0.0497(2) 


8 


0.0215(2) 


32 


0.2514(2) 


21 


0.1058(2) 


12 


0.0468(2) 


10 


0.0205(2) 


64 


0.2483(2) 


32 


0.1050(2) 


16 


0.0459(2) 


12 


0.0203(2) 


128 


0.2471(2) 


48 


0.1043(2) 


21 


0.0453(2) 


14 


0.0201(2) 


250 


0.2463(2) 


64 


0.1042(2) 


32 


0.0450(2) 


16 


0.0202(2) 


CO 


0.2458(2) 


CO 


0.1037(3) 


CO 


0.0446(4) 


CO 


0.0199(5) 



A site in the d-dimensional bcc lattice has 2 d neighbours, located at ^(±Xi ± 
X2 ± . . . ± Xd). This is dubbed the 'generalized bcc net' by Conway and Sloanei. 
A possible set of lattice vectors is b, = — ^(xj + . . . + Xj_i) + |(Xj + . . . + x^) for 
i = 1, . . . , d. The neighbours of a site are given in terms of these vectors by 

±b, 

±(b i -b j + b k ) 
±(bj - bj + b k - bi + b m ) 

The ci-dimensional diamond lattice is a lattice with a 2-point basis. Let us call 
the two points in the basis A and B. The full lattice is built by translation of the 
lattice basis over d independent vectors t j . Each A-site has d + 1 neighbours of 
type B. One of these neighbours is the B-site in the same basis, while the other 
d neighbours are the B-type sites at a relative location tj. Each B-site also has 
d + 1 neighbours, one of which is the A-site in the same basis and the others are 
the A-type sites at a relative location — tj. 

3. Lattice coding and cluster algorithm 

Two programs to calculate percolation thresholds were developed. One program 
was geared towards the handling of any desired lattice topology, the other towards 
speed and efficient memory usage for certain specific lattices' 1 . 

In the general purpose program to calculate percolation thresholds we used a 
generic method to specify a lattice. The bonds connected to lattice sites were 
coded in Fortran as NB_SITES(QMAX,NSITES) , where NSITES is the total number of 
sites in the lattice, and QMAX is the maximum number of bonds connected to a site 
(maximum coordination number). The integer value of NB_SITES(I ,S) was set to 
the bond number of the i th bond connected to site s. 

The sites connected to bonds were coded as NB_B0NDS (2 , NB0NDS) , where NB0NDS 
is the total number of bonds in the lattice. In other words, link i connects site 
NB_B0NDS(1,I) and NB_B0NDS (2,1). The two arrays NB_SITES and NB_B0NDS to- 
gether specify the topology of the lattice completely. 

a The programs, one in Fortran and one in C, can be obtained from the author. 



for i = 1, . . . , d, 

for i > j > k, 

for i>j>k>l>m, 



(3) 
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Table 3. The site percolation thresholds of d-dimensional diamond lattices as a function of linear 
lattice size L. The values in the last row are results of a fit to the scaling relation \p c (L)—p c (oo)\ oc 
L^ 1 /". Error estimates concerning the last digit are indicated between brackets. 





d = 2 




d = 3 




d = 4 




d = 5 




d = 6 
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Pes 


L 


Pes 


L 


Pes 


L 


Pes 


L 


Pes 


256 


0.6964(2) 


16 


0.4398(2) 


16 


0.3014(2) 


8 


0.2268(2) 


8 


0.1760(2) 


512 


0.6964(2) 


32 


0.4347(2) 


21 


0.2996(2) 


12 


0.2251(2) 


10 


0.1762(2) 


1024 


0.6969(2) 


64 


0.4317(2) 


32 


0.2989(2) 


16 


0.2248(2) 


12 


0.1770(2) 


2048 


0.6970(2) 


128 


0.4312(2) 


50 


0.2984(2) 


24 


0.2249(2) 


14 


0.1777(2) 


3072 


0.6970(1) 


250 


0.4306(2) 


64 


0.2983(2) 


30 


0.2251(2) 


16 


0.1783(2) 


oo 


0.6971(2) 


oo 


0.4301(2) 


oo 


0.2978(2) 


oo 


0.2252(3) 


oo 


0.1799(5) 



When we want to calculate a percolation threshold for such a lattice, we need a 
cluster algorithm, e.g. the one proposed by Hoshen and KopelmanLJ. This algorithm 
works by assigning a cluster label m to each site and bond. (It is assumed that 
to > 0.) In addition to the arrays to hold these labels, the algorithm uses one 
other array N(m). For each proper label, N(m) is greater than zero (by definition 
of 'proper'), in which case it holds the number of sites and bonds that belong to 
the cluster to. For a non-proper label JV(to) is negative, and —N(m) refers to a 
label to', which is the cluster that cluster to has been merged with. 

The algorithm consists of the following steps. 

• Initialise by setting N(m) <— for all to, and m next <— 1. 

• Loop over all sites and bonds of the lattice. 

• For each new occupied site and bond one encounters, one has to determine 
which cluster it belongs to. To this end, loop over the neighbours that one 
already determined a cluster label for. Three different situations can occur. 

=0 If there are no labeled neighbours, one has to define a new cluster label. 
Use set N(m) <— 1, and m n ext *— m ne xt + 1- 

=1 If there is one labeled cluster to, the current site/bond will get the label 
of that cluster. Set N(m) <- N(m) + 1. 

> 2 If there are two or more neighbours that belong to different clusters, then 
those clusters have to be 'merged': 

— Determine the lowest cluster label, mi ow , of these neighbours. 

— For all other labels m ther, set N(m ther) < nii ow , and 

N{m low ) <- N(m low ) + N(m other ). 

— For the current site or bond, set to <— to; oui , and N(m) <— N(m) + 1. 

In this algorithm, it is important to know the neighbours of a site or bond. The 
specification of neighbours is the only item that is lattice-specific. When one is 
interested in site-bond percolation, the neighbours of a site are directly given by 
the NB_SITES array, and the neighbours of a bond by the NB_B0NDS array. Therefore 
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the implementation of the Hoshcn-Kopclman algorithm is fairly straightforward for 
site-bond problems. There is but one important difference with the algorithm as 
described by Hoshen and Kopelman, which is that in their article, the networks 
could be traversed in an ordered way whereas here we cannot. Let us take as an 
example the square lattice in two dimensions. If the loop over sites is performed as 
DO I = 1, L 

DO J = 1 , L 

ENDDO 
ENDDO 

it is clear that neighbouring sites {i, j — 1} and {i — have been visited before 
site {i,j}, whereas {i + 1, j} and {i, j + 1} have not. 

Here, on the other hand, we are dealing with a generic lattice, which has no 
definite order, and hence one cannot traverse the lattice from one side to the other in 
a regular way. Therefore it is not a priori clear which neighbours have already been 
visited and hence have been assigned a cluster label. This problem can be solved 
by a suitable initialisation. For instance, we can define arrays to hold the cluster 
label per site (LS(NSITES)) and per bond (LB(NBONDS)). During initialisation we 
can set LS(I)=0 for occupied sites and LS(I)=-1 for empty ones, and similarly for 
bonds. In the couse of the algorithm, the sites for which LS(I)=0 will be assigned a 
valid label LS(I)=to. Since valid labels are greater than zero, one can, at any time 
during the algorithm, recognise which sites are 'empty', which ones are 'occupied 
but unassigned', and which ones are 'assigned'. 

Having performed the cluster algorithm, one still has to decide whether there 
is a percolating cluster or not. One has to define, therefore, what is percolation 
on this lattice: from where to where does a cluster have to extend, to be called 
percolating? In other words, one should define a beginning of the lattice and an 
end (or an 'IN' and an 'OUT'). For example, one can arrange the list of sites such 
that the last N out sites will be the exit of the lattice and the last N in before that 
will be the entrance. Having done so, one can loop over the entrance sites and mark 
the clusters that are connected to it. The final step is then to loop over the exit 
sites, to see whether there is a 'marked' cluster connected to the exit. If so, there 
is a percolating cluster, and otherwise there is not. 

The marking of a cluster can be done in various ways, for example by adding 
N s+ b to N(m), where N s+ b is the total number of bonds and sites. Because the 
total number of sites and bonds that can belong to any cluster is at maximum N s+ t,, 
this method is allowed: any cluster with N(m) > N s+ b is marked as connected to 
the entrance, and consists of N(m) — N s+ b sites and bonds. 

The cases of site-site percolation or bond-bond percolation are special cases of 
site-bond percolation. One can e.g. calculate the bond-bond percolation thresholds 
by simply making sure that all sites are 'occupied'. 

The special purpose program to calculate the site percolation thresholds for 
several lattices is both faster and more economic in the use of memory, and was 
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Table 4. The bond percolation thresholds of ci-dimensional Kagome, fee and bcc lattices as a 
function of linear lattice size L. The values in the row marked with oo are results of a fit to 
the scaling relation \p c (L) — p c (oo)| oc L -1 /" . Error estimates concerning the last digit are 
indicated between brackets. The values in the last row were obtained with an alternative method 
for extrapolation to L = oo (see text). 

d = 4 d = 5 
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Kagome 


fee 


bcc 


L 


Kagome 


fee 


bcc 


8 


0.1725(2) 


0.0473(2) 


0.0729(2) 


6 


0.1176(2) 


0.0227(2) 


0.0287(2) 


12 


0.1755(2) 


0.0484(2) 


0.0738(2) 


8 


0.1237(2) 


0.0237(2) 


0.0301(2) 


16 


0.1764(2) 


0.0487(2) 


0.0738(2) 


10 


0.1264(2) 


0.0252(2) 


0.0309(2) 


21 


0.1767(2) 


0.0490(2) 


0.0740(2) 


12 


0.1275(2) 


0.0255(2) 


0.0314(2) 


oo 


0.177(1) 


0.049(1) 


0.074(1) 


oo 


0.130(2) 


0.026(2) 


0.033(1) 


check 


0.179(2) 


0.049(1) 


0.075(1) 


check 


0.132(2) 


0.027(1) 


0.033(1) 



written in C. For the Bravais lattices, this is achieved by only using an array for 
the sites sites (x,y,. . .) and an array for N(m). These lattices can be traversed in 
a structured way, so that we know in advance which neighbours have been visited 
before, and which ones haven't. The specification of the neighbours of a site is done, 
for the four-dimensional fee lattice, by 

int ixp[] = { 1, 1, 0, 1, 0, 0, 0,0, 0, 1,1, 1>; 

int iyp[] = { -1, 0, 1, 0, 1, 0, -1, 0, 0,-1,-1, } ; 

int izp[] = { 0,-1,-1, 0, 0, 1, 0,-1, 0,-1, 0,-1 }; 

int iap[] = { 0, 0, 0,-1,-1,-1, 0, 0,-1, 0,-1,-1 }; 

These statements specify 12 neighbours, out of a total of 24 for a site on the four- 
dimensional fee lattice. The first six of these neighbours correspond with the first 
line in Eq. (Q), the last six with the second line. The above neighbours are only the 
ones that have been visited before, because the loop over sites is performed with 
the loop over 'ia' as the outer-most loop, and the loop over 'ix' as the inner-most 
loop. For each of the neighbours, ixp gives the x-displacement with respect to the 
current site. 

For non-Bravais lattices, like the diamond lattice, the sites array was given an 
extra dimension I, where / runs over the sites in the basis of the lattice. For the 
diamond lattice in d dimensions, the basis consists of 2 sites. The loop over I was 
implemented as the inner-most loop. The neighbours of a site can then be specified 
as follows. The sites of type A are treated differently from the sites of type B. 



int ilp[] 


= { 


1, 


1, 


1, 


1, 


-1 


} 


int ixp [] 


= { 


-1, 


o, 


o, 


o, 





} 


int iyp[] 


= { 


0, 


-1, 


o, 


o, 





} 


int izp[] 


= { 


0, 


o, 


-1, 


o, 





} 


int iap[] 


= { 


o, 


o, 


o, 


-1, 





} 



The first four neighbours are neighbours of the site of type A. They are sites of type 
B, located at a displacement of — x, respectively — y, . . . The last neighbour is a 
neighbour of the site of type B. It is the other site, type A, in the same basis. 
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4. Results 

The percolation thresholds of fee, bee, and diamond lattices are given in Tables @- 
||, for several lattice sizes. All the results were obtained with runs on a Sun Sparc 
workstation with 320 Mb internal memory. The values in Table ^ were calculated 
by use of the general purpose program. This program allowed only limited lattice 
sizes, for which it is not sure that the scaling relation 

\ Pc (L)- Pc (cx>)\ ocL- 1 /- (4) 

holds. Therefore an estimate was made of a possible systematic error. This was 
done by comparing a fit of the scaling relation to the last three data points and a fit 
to the last two data points. The estimate is given as the error margin for p c (oo) in 
the table. During the fitting procedure the value of the exponent v was fixed at 0.88 
in three dimensions, and v — 0.68, 0.57, 0.5 in 4, 5 and 6 dimensions respectivelyS. 

As a check on the results for p c (oo), we also used an alternative method to 
determine the value for the percolation thresholds at L = oo. When we define 



A(L) = VI < Pc{Lf > - < Pc (L) > 2 |, (5) 
we can use the relation! 



I < Pc(L) > — p c (oo)| oc A(L). (6) 

With the use of this relation, the value of p c (oo) can be fitted, without knowledge 
of the value of the exponent v. The results of this fit are shown in the row marked 
'check' in Table ||. They confirm the values of the row marked 'p c (oo)\ 

The results in Tables Q-|3| were obtained using the special purpose program, 
which was designed to suit the particular lattice at hand. Fits to the data points 
show that the values in the tables are within the scaling regime, see Fig. |l|. Therefore 
the error margins quoted here for p c (oo) are the ones indicated by the fit. The values 
for the fee lattices in Table |l| are different from those reported by Zallero (0.098 
in d — 4 and 0.054 in d = 5). However, the cross-check with the general purpose 
program confirmed the numbers of Tables 0-^. Furthermore, it is not clear how the 
Zallen results have been calculated, nor how large their estimated error margin is. 



5. Discussion 

Percolation thresholds of many lattices are listed in Table ||. The results for the 
deviation of C = "Pcb lrom 1 i s a l so listed. When we focus on the results 

for the invariant C, it is interesting to compare the triangular lattice with the 
octagonal lattice. Both lattices are isotropic lattices, with (average) coordination 
number 6. Roth lattices are fully triangulated, and hence have equal site percolation 
thresholdsO p C8 = -|. Nevertheless, their bond percolation threshold differ: Ueb = 
2sin(vr/18) = 0.347 296 . . .1 for the triangular lattice vs. p cb = 0.3237 ±O.OOO60for 
the octagonal lattice. Therefore, the value of C for these lattices is quite different, 




10 3 10 4 10 5 10 6 10 7 

lattice size N = L 

Fig. 1. The scaling of the bcc site percolation thresholds as a function of the lattice size 
see Eq. (^|). The fits to the data have a slightly different slope, because the exponent v is 
different for each dimension up to d — 6. 



namely C — 0.97 for the triangular and C = 0.88 for the octagonal lattice. Also 
for the dual of the octagonal lattice, the value for C deviates substantially from 
unity: C = 0.91. 

One can also compare the cubic lattice and the 3-dimensional Kagome la 
tice. The latter has a much higher site percolation threshold (0.3895 ± 0.0002L 
vs. 0.311 604 ± 0.000 00060), but the bond percolation thresholds lie closer to- 
gether (0.2706 ± 0.00090 vs. 0.248 812 6 ± 0.000 000 50). Although both lattices 
are isotropic and have coordination number 6, the respective values for C are 0.99 
and 0.75. 

The dual of the diamond lattice and the 3-dimensional Kagome lattice form 
a special pair too. It can be shown that the site percolation thresholds of these 
lattices are equal, although the lattices have a different coordination number!^. On 
the other hand, their bond percolation thresholds are different, which is reflected 
in the values for C being 0.75 and 0.65. Note further that the dual of the bcc 
lattice, which was not incorporated in the Galam and Mauger studyB, also shows a 
substantial deviation, with C — 0.89. 

Moreover, in high dimensions the dependence of percolation thresholds on the di- 
mension d is not universal. For the proposed invariant, Galam and Mauger assumed 
that all thresholds have equal scaling behaviour, such as the relation p c oc l/(2d— 1) 
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site percolation 

o 
□ 



0.01 



10 



100 



(d-1)*(q-1) 

Fig. 2. Site percolation thresholds as a function of (d — l)(q — 1), see Eq. ([?]) 



that was established some time ago for hypercubic latticesO. However, it has been 
shown recently!^ that there are d- dimensional Kagome lattices with a different scal- 
ing behaviour, namely p c oc 1/d. For these lattices the values for C decrease with 
dimension, until C — 0.58 for d = 5. Also for the fee, bec, and diamond lattices, 
the behaviour as a function of dimension has not yet been captured fully. For each 
of the lattices the C-value deviates more strongly from unity, as the dimensional- 
ity increases. The scaling of nercolation thresholds as a function of dimension was 
fitted by Galam and Maugertl to 



p„oc{(d-l)(g-l)}- 



(7) 



for site percolation thresholds, and 



Pcb oc 



(d- 



(8) 



for bond percolation thresholds. The thresholds have been plotted in Figs, g and 3, 
together with the Galam and Mauger fit. The bond percolation thresholds, Fig. 3, 
all seem to have the same trend, and the fit captures that trend quite well. The site 
percolation thresholds, however, show a much more complicated behaviour (Fig. |^). 
The deviations from the Galam and Mauger fit are increasing with dimension, but 
also it looks as if there is a spread in the thresholds, even though they follow a 
general trend. The Kagome lattices should be considered separately, since thei 
scaling behaviour has been shown to differ from the scaling for hypercubic latticed 



I 
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0.01 



10 

(d-1)*(q-1)/d 



Fig. 3. Bond percolation thresholds as a function of (d — l)(q — l)/d, see Eq. 



In conclusion, site and bond percolation threshold have been calculated for sev- 
eral lattices in up to six dimensions. The scaling behaviour of these percolation 
thresholds as a function of dimension is sometimes different than for other known 
lattices. It may therefore be necessary to introduce more separate classes to al- 
low for the variation in behaviour, e.g. a separate class for each dimension. If 
the deviations would then become smaller again, the invariant proposed by Galam 
and Mauger would still have the remarkable property of being independent of the 
coordination number. 
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Table 5. Site (pes) and bond (p c b) percolation thresholds for various lattices and their dual lattices. 
The deviation of the quantity C = Sp c J~^ a ' p 1 J b ab /d from 1 is listed under C — 1. Between brackets 
are error estimates concerning the last digit. The thresholds without error margin are <J*fte% 
For the site percolation thresholds on cubic lattices in higher dimensions, see also Refs. EJ'Ej. 
The lattices in the upper part of the table are two-dimensional, including a number of aperiodic 
lattices; the ones in the middle are three-dimensional, and in the lower part of the table the higher 
dimensional lattices are listed. 



lattice 



Pes 



Pcb 



C — 1 pes dual 



p ch dual 



C-l 



square 
Kagome 
pentagonal 
bowtie 
triangular 
octagonal 

Penrose 
octag-chem 

octag-ferr 
dodec-chem 

dodec-ferr 

cubic 
3d-Kagomc 
diamond 
trian. stack 
bec 
fee 
hep 

4d-cubic 
4d-Kagome 
4d-diamond 
fee ( l ZV) 
4d-bcc 

5d-cubic 
5d-Kagomc 
5d-diamond 
fee ('D5') 
5d-bcc 

6d-cubic 
6d-Kagome 
6d-diamond 
fee ('IV) 
6d-bcc 



0.592 746 OH 1 
0.652 703 
0.6471(6) 
0.5 




0.5837(3)0 
0.585(1)1 
0.543(2)55 
0.628(2)22 
0.617(3)E3 
0.311 604(6)liil 
0.3895(2)0 
0.4301(2)r-| 
0.2623(2) }£ 
0.2458(2) UJ 
0.1994(2) \£ 
0.1990(2)E 

0.196901 fe )B 
0.2715(3)0 
0.2978(2) 
0.0842(3) 
0.1037(3) 

0.1407(3) \ 
0.2084(4)E 
0.2252(3) 
0.0431(3) 
0.0446(4) 

0.1079(5)1 
0.1677(7)E 
0.1799(5) 
0.0252(5) 
0.0199(5) 



0.5 13 n 

0.524 405p^3)EIl 
0.5800(6)li3 n 
0.404 518... p 
0.347 296it 
0.3237(6)i 

0.4770(rfr; 
0.478(3) U 
0.402(5) y 
0.538(1)0-1 
0.4950(5)B 

0.248 812i#<5)y 
0.2709(6)^ 
0.3893(2)^ 
0.1859(2)0 n 
0.180 287 5(10) yj 
0.120 163iS^10)E3 
0.1199(2)i 

0.1600(2)E 
0.177(1) r, 
0.2715(3)0 
0.049(1) 
0.074(1) 

0.1181(2)H 
0.130(2) i-i 
0.2084(4)0 
0.026(2) 
0.033(1) 
0.0943(2)H 



0.1677(7 



0.04 
-0.02 

0.01 
-0.05 
-0.03 
-0.12 

0.01 
0.01 
-0.04 
-0.01 
0.01 

-0.01 
-0.25 
-0.06 
-0.03 
0.04 
-0.05 
-0.05 

-0.01 
-0.35 
-0.12 
0.13 
0.25 

-0.04 
-0.42 
-0.16 
0.39 
0.69 

0.01 
-0.19 



0.5848(2) tS 

0.5502(8) 

0.6653(6) 

0.6971(2) 

0.7297(4) 

0.6381(3) 



0.3895(2) 
0.3701(2) 
0.4560(6) 
0.3341(5) 
0.3101(5) 



0.475 594^3)0 
0.4196(6)E3i-i 

0.595 481.. 
0.652 703r 
0.6771(6) 

0.5233(2) 



0.2350(5) 
0.3093(2) 
0.4031(6) 
0.2703(3) 
0.2573(3) 



0.00 
-0.01 
-0.02 
-0.02 
-0.09 

0.01 



-0.35 
-0.05 
-0.11 
-0.03 
0.04 



